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Abstract 

For high precision measurements of K decays, the presence of radiated photons 
cannot be neglected. The Monte Carlo simulations must include the radiative cor- 
rections in order to compute the correct event counting and efficiency calculations. 
In this paper we briefly describe a method for simulating such decays. 
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1 Introduction 



Many measurements on K decays have reached a statistical error close to 1% 
or less. With such precision, the presence of radiated photons, and in general 
the effect of radiative corrections, cannot be neglected. Furthermore, the 
treatment of the radiative corrections is explicitly required for the extraction 
of many physical quantities, such as the CKM matrix element V us and the 
phase shifts 5q — 62, at precisions of a few percent (or above). Hence, it is 
mandatory to include the effect of radiated photons in the Monte Carlo (MC) 
simulations. 
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Figure 1: Energy spectrum for K° — > irewy MC decays. 2.4% of the events 
have £L > 30 MeV. 



There are two aspects of a measurement that are affected by the pres- 
ence of radiated photons [T|: the geometrical acceptance and the counting 
of the events. About 2.4% of K° — > irew-f decays have a photon with an 
energy above 30 MeV, as shown in Fig. Q corresponding to about 10% of 
the center-of-mass energy. These photons soften the momentum spectra of 
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the electron and pion that are actually detected in an experiment, changing 
the geometrical acceptance and the distributions of kinematic quantities of- 
ten used to select or count the number of signal events. If these effects are 
neglected, errors up to few percent can result in the measurement of a given 
branching ratio j2] [3]. 



2 Bremsstrahlung and infrared divergences 

The main problem in simulating radiative decays is the presence of infrared 
divergences: the total decay width for single photon emission, computed at 
any fixed order in a, is infinite. A finite value is obtained only by summing 
the decay widths for the real and virtual processes calculated to the same 
order in a. As shown in Ref. in the limit of soft photon energy, we can 
"re-sum" the probabilities for multiple photon emission to all order in a. The 
rate for the decay process i — > / accompanied by any number of soft photons 
with total energy less than J5 7 is given by: 

r ind (E 7 ) = r (^)V + o{b 2 ) + o(£ 7 )). (i) 

Here T is the unphysical decay width for the process i — > / without final 
state photons, and b is a function of the particle momenta, is positive and of 
order a; it is given by: 

b = ~ VmVne m e n (3 mn In 1 * ^ mn , (2) 

m,n 1 Pmn 

where m and n run over all the external particles, e n is the charge of the 
particle n, rj = +1 or — 1 for an outgoing or incoming particle, and (3 mn is 
the relative velocity of the particles n and m in the rest frame of either: 
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A is an energy cut-off that can be chosen as the mass M of the decaying 
particle. 

Differentiating rj nc /(i£ 7 ) with respect to _E 7 , we obtain an integrable dif- 
ferential distribution: 

dE~, ~ M b ~ dEy \m) ' [ ' 

where we have neglected second order terms in b. To 0(a) this can be 
identified with the single-photon emission probability, and indeed, it can be 
written in terms of: 

Brem p b 

dE~f E-y 

The presence of the extra factor (E 1 /M) b ensures the integrability of Eq. (j3J) 
in the limit E 1 — > 0. 

In the derivation of Eq. (0), in Ref. no explicit integration is required 
on the momenta of particles other than those of the photons. The result in 
Eq. (jlj) can thus be applied to differential decay withs: 



d^ incl dY Brem ( E-y \ ^ 

dE~d~l ~ ~dE~dJ_ KM J 



(6) 



where £ represents the independent kinematic variables of the decay process 
without photons, and where: 

dT B rem _ d£o K£) 

dE^d£ d£ E 7 ' U 

Note that, while for two body decays b is a constant, for decays with more 
particles in the final state, the velocities j3 in Eq. Q and thus b, depend on 
the variables £. 



3 MC simulation for radiative decays 

While for a complete MC simulation we need the decay width for all values 
of Ej, the relation in Eq. (jHJ) is true for only soft-photon emission. However, 
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the value of the exponent b, which is of order a, is about 0.01. Hence, 
while (E^/M) b is a big correction for E 1 — > 0, its value is close to 1 when 
E 1 — ► M. Therefore, if we use the complete differential decay width at order 
a for the emission of a photon, dT Breml 'dE^d^, instead of its approximation 
for low energies in Eq. (J7J), the decay width in Eq. (0) represents a good 
approximation for the entire energy spectrum. 

In summary, the "recipe" for writing the amplitude for the process % — > f'y 
is the following: 

1) Calculate the b factor using Eq. (0). 

2) Calculate the amplitude Mj^/ 7 at order a. 

3) Fix the divergence in the squared amplitude by multiplying by (E^/M) b : 

\^-> h \l oIR =\M^ h \ 2 {E,/M) b . (8) 

3.1 MC generators 

Following the recipe outlined in the previous section we have written the MC 
generators for the decays: K° — > 7r7r7, K° — > irei"-/, K° — > TT^w-f, K° — > 
— > 7T7T7, — > ez/7, i ; T ± — * /i^7, -ft' 1 ' 1 — > rrewy, — > TV/jLVJ, 
and — ► We obtained the amplitudes M^f^ at order a mainly 

from Ref. jH], where they are calculated using chiral perturbation theory 
at order p 2 and p 4 . Concerning the semileptonic modes, for simplicity we 
have used only the p 2 expression of the Mj_»j 7 amplitudes. At this order, 
the form factors /+ and /_ are equal to 1 and respectively. In order to 
take into account the leading dependence on the variable t = (px — P 7 ^) 2 we 
have multiplied the amplitude by an overall factor (1 + X + t/M 2 + ). 1 The 
uncertainty related to this approximation is discussed in the following. 

Large MC production in high energy experiments puts stringent limits on 
the time needed to generate one single event. This time should not exceed 
: We used the value A + = 0.03. 
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the time needed to track the particles inside the detector. In the KLOE 
experiment this time is about a few milliseconds. We used a combination of 
MC sampling techniques, the acceptance-rejection (Von Neumann) method 
and the inverse transform method 0, to reach this goal JU]- The average 
time for generating one event is a fraction of a millisecond. 

We compared the fraction of events with a photon above an energy thresh- 
old predicted by the MC simulation, with theoretical expectations and ex- 
perimental results, for several decay channels of the kaon. For instance, the 
MC prediction for the fraction of K° — > 7r + 7r~7 in which a photon has energy 
greater than 20 MeV (50 MeV) is equal to 7.00 x 10~ 3 (2.54 x 10~ 3 ), in agree- 
ment with the measured value (7.10 ± 0.22) x 10~ 3 ((2.56 ± 0.09) x 10~ 3 ), 
and that from theoretical predictions 7.01 x 10 -3 (2.56 x 10~ 3 ) for the Kg 2 
(see Ref. 0). 

Moreover, we calculated the ratio 

D r(i^ e37 , E 1 > E, 6 n > 6) 

R = ip^w) ' (9) 

giving the fraction of decays k —>■ nevy with a photon with energy E 1 above 
E and angle between the electron and photon 9 ei above 9, for different values 
of the energy E and angle 9, and we compared it with theoretical predictions 
and, when ever possible, with experimental values. MC calculations of R° 
for K° — > newy decays are shown in Tab. ^ while theoretical expectations 
from Ref. ^T] and Ref. |X2] 3 are shown in the top and middle part of Tab. El 
respectively. Since the decay width in the denominator of Eq. Q is inclusive 
of the photon emission, following Ref. ^5], we divided the predictions of 
Ref. [TT] and [12] by a factor (1 + 51), where 5% is the total electromagnetic 
correction extracted from Ref. ^3] • Note that a different approach is used in 
Ref. [7j, where most of the electromagnetic corrections are absorbed in /+(0) 
and therefore cancel out in the ratio R°. 



2 These numbers refers only to the inner bremsstrahlung (IB) term. 

3 We quote the values from Ref. JE], obtained for A+ = 0.03, as in our simulation. 
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The results for R° from a recent MC simulation, described in Ref. [T4] . 
are shown in the lower part of Tab.|21 Experimental results for R° have been 
recently published in Ref. [T^j by the KTeV Collaboration, for two values of 
the photon energy and angle: 

T(K e3 ~, > 10 MeV, 6^ > 0°) , 

' n - '- = (4.942 ± 0.062) x 10~ 2 , 



r(^ e3 (7)) 

ft -> 0C\°\ 

= (0.916 ± 0.017) x 10" 



T(K e37 , E, > 30 MeV, 9 n > 20°) , n ft1 « , n ni ~ M , n _ 2 



r(^ e3 (7)) 

and in Ref. [TSj by the NA48 Collaboration: 

r(iT e 3 n ,, EL > 30 MeV, 9 e -* > 20°) . , nmiN , 

1 7 \ (K ( » 7 " = (0.964 ± 0.008^) xlO- 2 . 

r(A e3 (7)) 

MC calculations and theoretical expectations of R^ for — > 7rez/7 
decays are shown in Tab. El and in Tab. 0] Since the radiative correction 8% 
for T(K~f 3 ) is compatible with zero, 8% = +0.06(20)%, it has been neglected. 

In Tab. ^ and Tab. E] the first error is statistical while the second is 
the systematic one. We computed the systematic error by comparing the 
difference between the branching ratios calculated at order 0(p 2 ) and at 
order 0(p A ) in Ref. jHj, with the variations of R we observed in the simulation 
when the term (1 + A+t/M 2 + ) is included, or not, in the decay amplitude. 
In Ref. jH] the branching ratios evaluated for E^ > 30 MeV and # 7 > 20° 
increase by about 6% in the 0(p 4 ) calculation, while the variations of R in 
the simulation are about 3%. Hence, we used the variations observed in the 
simulation for each of E~ and # 7 as systematic errors. Moreover, we changed 
the value of the cut-off A from Mr to Mr-/2 to check the stability of the 
results, resulting in negligible variations of the order of the statistical errors 
shown in the two tables. 

For both R° and i? ± the absolute differences between the MC simulation 
and theoretical predictions are below 5 x 10 -4 , and below the quoted sys- 
tematic errors. Moreover, such errors are smaller than the relative errors on 
measured branching ratios of kaon decays j2] [2] 0j. 
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#/E(MeV) 


10 


20 


30 


40 


0° 


4.908(7)(56) 


3.252(6)(47) 


2.364(5)(43) 


1.784(4)(38) 


10° 


2.450(5)(38) 


1.657(4) (36) 


1.223(3)(31) 


0.937(3)(28) 


20° 


1.864(4) (33) 


1.262(4)(30) 


0.933(3)(28) 


0.716(3)(24) 


30° 


1.516(4)(30) 


1.028(3)(28) 


0.760(3)(25) 


0.583(2)(22) 


40° 


1.264(4)(27) 


0.859(3)(25) 


0.636(3)(23) 


0.488(2)(21) 


50° 


1.066(3)(24) 


0.726(3)(23) 


0.538(2)(20) 


0.414(2)(19) 



Table 1: Ratios R° x 10 2 for K°e3 decays obtained in the MC simulation 
with 10 7 events. The first error is statistical while the second is systematic. 



Ref. [II. 9/E(MeV) 


10 


20 


30 


40 


0° 


4.94 


3.25 


2.35 


1.77 


10° 


2.48 


1.67 


1.24 


0.95 


20° 


1.90 


1.29 


0.95 


0.73 


30° 


1.55 


1.05 


0.78 


0.60 


40° 


1.29 


0.85 


0.66 


0.50 


50° 


1.10 


0.75 


0.56 


0.43 


Ref. 12j e/E(MeV) 


10 


20 


30 


40 


0° 




3.27 


2.37 


1.78 


Ref. 14J 6/E(MeV) 


10 


20 


30 


40 


0° 


4.93(6) 




2.36(3) 




20° 


1.89(2) 




0.96(1) 





Table 2: Ratios R° x 10 2 for K°e3 decays: (Top) Values listed in Ref. jH] ; 
(Middle) Values obtained from Ref. ^2] ; (Bottom) Values from the MC sim- 
ulation in Ref. [T3]. The values from Ref. ^T] and Ref. ^2] have been mul- 
tiplied by (1 + SI)' 1 where 5% = 1.04(20)% from Ref. |T3]. 
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#/E(MeV) 


10 


20 


30 


40 


0° 


4.223(6)(34) 


2.825(5)(28) 


2.069(5)(27) 


1.572(4)(25) 


10° 


1.781(4)(18) 


1.238(3)(17) 


0.936(3)(17) 


0.732(3)(15) 


20° 


1.211(3)(16) 


0.854(3)(15) 


0.652(3)(14) 


0.515(2)(13) 


30° 


0.881(3)(15) 


0.630(3)(14) 


0.488(2)(13) 


0.388(2)(12) 


40° 


0.659(3)(14) 


0.477(2)(13) 


0.373(2)(11) 


0.300(2)(10) 


50° 


0.493(2)(12) 


0.363(2)(11) 


0.289(2)(11) 


0.234(2)(10) 



Table 3: Ratios R^ 1 x 10 2 for K e3 decays as obtained in the MC simulation 
with 10 7 events. The first error is statistical while the second is systematic. 



Ref. [n #/E(MeV) 


10 


20 


30 


40 


0° 


4.24 


2.80 


2.03 


1.53 


10° 


1.78 


1.22 


0.91 


0.71 


20° 


1.20 


0.83 


0.63 


0.50 


30° 


0.87 


0.61 


0.47 


0.37 


40° 


0.65 


0.46 


0.36 


0.29 


50° 


0.48 


0.35 


0.27 


0.22 


Ref. 12 6/E(MeV) 


10 


20 


30 


40 


0° 




2.82 


2.04 


1.54 



Table 4: Ratios R± x 10 2 for K^eS decays: (Top) Values listed in Ref. P]; 
(Bottom) Values obtained from Ref. [12]. The radiative correction 5% = 
+0.06(20)% is compatible with zero and has been neglected. 
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4 Conclusions 



We overcome the problem of infinite probabilities in radiative processes by 
extending the soft-photon approximation of Ref jHj to the whole energy range. 
The spectra produced with MC generators developed with this technique, 
agree well with other theoretical calculations and with available experimental 
data. The systematic error could be further reduced by using full 0{p A ) 
calculation for the amplitudes, or even adding the 0(p 6 ) results recently 
published in Ref. [7j. 

We used MC sampling techniques to reduce the time needed to generate 
one event below 1 ms. The MC generators, routines written in Fortran, 
have been included in the official KLOE library and have been used for 
large MC productions. Such routines are available on request to the author 
(claudio.gatti@lnf.infn.it) . 
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